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Abstract .  A  case  study  is  presented  of  a  hydrologic  investigation 
the  Red  Run  Drain-Lower  Clinton  River  watershed,  an  area  near 
Detroit,  Michigan,  that  has  undergone  urbanization  since  the 
1940's.  The  purpose  of  the  study  was  to  determine  peak-discharge 
frequencies  at  gaged  and  ungaged  locations  for  existing  and 
future  conditions.  Population  density  was  used  as  an  indicator  of 
urbanization  in  relationships  defining  unit  hydrograph  parameters 
and  hydrologica 11 y  significant  impervious  area.  Input  parameters 
for  a  single  event  rainfall-runoff  simulation  model  (HEC-1)  wer<’ 
developed  to  reflect  watershed  conditions  in  the  years  1940, 
1950,  1960  and  1975.  The  input  parameters  were  verified  by 

reconstructing  observed  flood  events  that  occurred  at  these 
points  in  time.  Sets  of  synthetic  winter  and  summer  storm 
hyetographs  were  input  to  HEC-1  to  develop  a  series  of  curves  for 
two  gaging  stations  that  relate  peak  discharge  to  magnitude  of 
synthetic  storm  for  each  watershed  condition.  The  curves  were 
used  to  transform  the  series  of  recorded  annual  peak  discharges  at 
each  gage  to  a  stationary  series  that  reflects  1975  watershed 
conditions.  Discharge  frequency  estimates  were  then  developed 
for  ungaged  locations  using  winter  and  summer  synthetic  storms 
that  were  assigned  exceedance  frequencies  consistent  with  actual 
exceedance  frequencies  at  the  gaged  locations.  Projections  of 
future  population  density  were  the  basis  for  developing  HEC-1 
input  parameters  representing  year  2000  and  2025  watershed  con¬ 
ditions.  Estimates  of  peak  discharge-frequencies  for  the  future 
conditions  were  made  at  the  gaged  and  ungaged  locations  using  the 
methods  described  above. 


Introduction 

Determining  the  impact  of  urbanization  on  flood 
flow  frequency  is  an  area  of  substantial  interest  for 
flood  control  planning  and  design  and  for  studies 
required  by  the  national  flood  insurance  program. 
This  paper  reports  methodology  used  in  a  Corps  of 
Engineers  study  (7)  to  evaluate  channel  modification 
proposals  for  the  Red  Run  Drain,  which  is  located 
immediately  north  of  Detroit  in  the  Clinton  River 
basin.  The  study  was  conducted  to  determine  peak- 
discharge  frequency  curves  at  two  gaged  and  numerous 
ungagt?d  locations  for  existing  (1975)  land  use  con¬ 
ditions  as  well  as  for  conditions  in  the  years  2000 
and  2025.  Methodology  is  reported  herein  for  con¬ 
verting  a  nonstationary  time  series  of  annual  peak 


discharges  to  a  stationary  time  series  representing  a 
specific  set  of  land  use  conditions,  and  procedure, 
are  described  for  developing  peak-discharge  frequency 
curves  at  ungaged  sites. 

Previous  Studies 

A  number  of  studies  relating  to  the  hydrology  of 
South*' astern  Michigan  have  been  reported  (1,  2,  3,  4, 
5,  6).  In  particular,  Brater  and  Sherrill  (2)  pro¬ 
vided  the  results  of  an  extensive  project  initiated 
in  1964  to  develop  a  procedure  for  estimating 
frequency  of  storm  runoff  for  small  basins  in  various 
stages  of  urbanization .  Results  of  that  project  wer* 
used  in  a  number  of  instances  in  the  study  desci'b.'d 
in  this  paper. 


♦Paper  presented  at  the  International  Symposium  on  Urban  Storm  Runoff 
University  of  Kentucky*  Lexington,  Kentucky,  July  23-26  1979. 
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Watershed  Characteristics  rection  for  11  miles,  through  the  City  of  Mt . 

'  Clemens,  to  its  outlet  in  Lake  St.  Clair. 

The  Clinton  River  Watershed  covers  an  area  of 

760  square  miles  in  the  southeastern  portion  of  the  Four  miles  upstream  of  the  Forks  Area  the  Main 

lover  peninsula  of  Michigan.  (Figure  1)  The  Clinton  Branch  is  joined  hy  the  Red  Run  Drain.  The  Red  Run, 

River  Watershed  includes  both  rural  and  urban  areas,  which  drains  an  area  of  approximately  140  square 

and  roughly  resembles  a  maple  leaf,  with  the  stem  or  mites,  originates  in  the  heavily  urbanized  northern 

base  the  outlet  to  Lake  St .  Clair.  The  Main  Branch  of  suburbs  of  Detroit,  and  flows  in  a  northeasterly 

the  Clinton  River  originates  in  the  hilly  country  of  direction  through  the  rapidly  urbanizing  sections  of 

Oakland  County,  vhert  a  great  number  of  lakes,  rang-  western  Macomb  County.  (Figure  2) 

ing  in  size  from  a  few  acres  to  several  square  miles, 

remain  a a  a  legacy  of  the  last  period  of  glaciation.  Three  distinct  types  of  topography  are  evident 

It  is  from  these  lakes  that  the  river  obtains  much  of  the  Clinton  River  Watershed:  The  hilly,  lake 

its  base  flow.  The  Main  Branch  flows  in  a  general  covered  region  that  comprises  the  head  waters  of  the 

southeasterly  direction  for  69  miles,  through  several  Main  Branch;  the  rolling,  well-drained  country  that 

sam  Her  urban  areas  (Pontiac,  Rochester,  Utica),  comprises  the  watershed  of  the  Middle  and  Worth 

where  it  is  joined  by  the  North  and  Middle  Branches,  Branches;  and  the  flat,  poorly  drained  areas  that 

which  drain  relatively  rural  areas  of  201  and  77  extend  in  a  broad  plain  along  Lake  St.  Clair  and 

square  miles,  respectively.  The  confluence  of  the  comprise  the  watershed  of  the  Red  Run.  ft  is  in  this 

three  tributaries,  known  locally  as  the  "Forks  Area,"  area  of  poorly  drained  glacial  rill  that  the  majority 

produces  the  Clinton  River  proper,  which  meanders  development  in  the  watershed  has  taken  place.  The 

through  a  level  plain  in  a  generally  easterly  di-  Drain  and  it  tributaries  have  been  ex- 


Flgure  1.  Clinton  River  Watershed 


tensive ly  modified  to  provide  adequate  capacity  for  aval l  able  for  tin*  period  1960-19/3  from  the  South* 

•torn  drainage  from  the  urban  areas.  However,  the  eastern  Michigan  Ra ingage  Network,  of  which  about  25 

rate  of  urbanisation  has  exceeded  the  rate  of  stations  lie  in  and  around  the  Clinton  River  Water- 

drainage  improvement  in  recent  years  and  wide-spread  shed.  Streamflow  records  available  for  the  30  U.S. 

overland  and  basement  flooding  has  resulted.  Geological  Survey  gaging  stations  on  the  Clinton 

River  and  tributaries  extend  to  a  maximum  period  of 
Su»fim  /  of  Available  Data  41  years  for  the  Main  Branch  Clinton  River  gage  at 

Moravian  Drive  in  Mt.  Clemens. 

Meteorologic  and  streamflow  records  for  the 
Clinton  River  Basin  are  relatively  comprehensive  and 

provide  excellent  coverage  of  the  more  recent  floods.  Parameter  Development  and  HfiC-1  Calibration 

Precipitation  records  for  the  City  of  Mt*  Clemens, 

considered  typical  of  the  basin,  are  available  from  Basin  modeling  consists  of  describing  the  topo- 

about  1897  to  the  present.  Hourly  rainfall  data  are  logic  structure  of  a  basin  (drainage  basin  bound- 


LEGEND: 

•  SEMCOG  RAIN  reO  RUN  DRAIN-LOWER  CLINTON  RIVER 

GAGING  STATIONS 
_  U.S.6.S.  STREAM 

GAGING  STATION  RED  RUN  BASIN 

H?)  SUB  BASIN  NUMBER 
sJ&Vs  WATER  COURSE 


SUB  BASIN  80UN0RY 
FREQUENCY  CURVE 


Table  1.  Listing  of  U.S.G.S.  Gages 


CAGE 

STATION 

TEARS 

OF 

RECORD 

04161000 

Mato  Branch 

Clinton  R.  4  Auburn  Hta. 

19 

04161100 

Galloway  Cr-  nr.  Auburn  Hta. 

16 

04161340 

Paint  Cr.  at  Rochcatar 

16 

04166000 

Stoney  Cr.  nr.  Uashlnjton 

17 

04162900 

Rad  Run  Watarahad 

Blf  Baavar  nr.  Uarren 

17 

0416)400 

Plus  Brook  at  Utica 

10 

04164000 

Clinton  Rivar 

Clinton  Rivar  nr.  Fraser 

20 

04164500 

N.  Br.  Clinton  Rivar  nr.  Mt.  Clenana 

28 

04163500 

Clinton  Rivar  at  Ht.  Clenana 

41 

aries,  stream  channels,  and  the  logical  relationships 
between  the  drainage  areas  and  stream  channels)  and 
defining  the  parameters  that  describe  the  precipita¬ 
tion-runoff  response  of  the  subbasins  and  the 
channels  that  comprise  the  basin.  The  basic  driving 


input  to  the  model  is  a  time  history  of  precipitation 
for  each  subbasin  for  the  duration  of  the  storm  event 
mder  study.  The  basic  output  is  computed  flood 
hydrographs  that  represent  a  time  history  of  the 
hydrologic  response  of  the  basins  and  stream  channels 
under  study. 

Computer  Program  HEC-1  (8)  uses  a  lumped-para¬ 
meter  approach  for  determining  runoff  from  each  sub- 
area  under  cons iderat ion,  i.e*,  the  parameters  and/or 
computations  for  a  particular  subbasin  are  assumed  to 
apply  for  its  entire  area.  Because  of  the  many  points 
of  interest  throughout  the  Red  Run  basin  and  the 
Lower  Clinton  River,  the  study  area  was  divided  into 
a  total  of  22  subbasins  ranging  in  size  from  roughly  5 
to  200  square  miles  as  shown  in  the  schematic  diagram 
in  Figure  3. 


Unit  Hydro^raph  Parameters  -  To  develop  the  required 
unit  hvdrographs  for  past,  existing  and  future  land 
use  conditions  for  the  numerous  subbasins  of  the 
Clinton  River  Watershed,  regional  relationships  were 
determined  for  unit  hydrograph  parameters.  In  ad¬ 
dition,  unit  hydrograph  parameters  for  five  rural 
subbasins  were  determined  by  optimization  techniques 
From  recorded  data. 


Figure  i.  II Sr . . Hie  Diagram 
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Table  2,  Drainage  Basin  Characteristics  and  Unit  Hvdr^grnph  Parameters 


Drainage 

Basin 

esc:; 

Cage 

No. 

Ar<  a 

Sn.  Mi. 

i  *1  h<li 

1)4  if,  i  1  y 
IVrM>ii.s/ 

So.  Mi. 

<vu) 

CFS 

.r  (2) 
fK 

Hrs 

Big  Beaver  Creek  near 
Warren.  MI 

4-1629.0 

23.5 

774 

799 

13.0 

Plum  Brook  near 

Utica,  MI 

4-1633.0 

22.9 

731 

779 

15.9 

Rouge  River 

Birmingham,  MI 

4-1660.0 

36.9 

1212 

1 107 

12.5 

Rou^e  River 

Southfield,  MI 

4-1661.0 

88.0 

1390 

2288 

15.0 

Rouge  River 

Detroit,  MI 

4-1665.0 

185.0 

2710 

3700 

27-0 

Clinton  River 

Mt.  Clemens,  MI 

4-1655.0 

734.0 

1000 

9836 

38.0 

McBride  Drain  near 
Macomb,  MX 

4-1644.5 

5.8 

114 

223 

9.3 

Tupper  Brook 

Ray  Center,  MI 

4-1642.5 

8.6 

60 

318 

5.0 

Evans  Ditch 

Southfield,  MI 

4-1662.0 

9.5 

2376 

603 

17.0 

East  Branch  Coon  Cr. 
Armada,  MI 

4-1643.0 

13.0 

56 

520 

10. 0 

Cloede  Ditch  near 
Waldenburg,  MI 

4-1652.0 

16.0 

560 

432 

14.1 

Upper  Rouge  River 
Farmington,  MI 

4-1663.0 

17.5 

1129 

700 

10.0 

Red  Run  near 

Royal  Oak,  MI 

36.5 

7500 

5110 

4.0 

(1)  Qp  =  unit  hydrograph  peak  discharge. 

(?.)  Tg  ■  period  of  rise  of  unit  hvdrograph. 


Criteria  for  relating  unit  hydrograph  parameters 
to  watershed  characteristics  for  use  in  southeastern 
Michigan  were  given  by  Brater  and  Sherrill  (2).  Unit 
hydrograph  parameters  for  13  gaged  basins  in  south¬ 
eastern  Michigan  are  summarized  in  Table  2.  The  unit 
hydrograph  parameters  in  Table  2  are  for  average  unit 
hydrographs  developed  from  a  number  of  storms  for 
each  basin.  Extensive  analysis  by  Brater  and  Sher¬ 
rill  showed  that  population  density  and  drainage  area 
were  the  most  significant  watershed  characteristics 
in  determining  unit  hydrograph  parameters.  Data  from 
Table  2  were  entered  into  HEC's  multiple  Linear  Re¬ 
gression  computer  program  to  develop  equations  for 
unit  hydrograph  peak  discharge,  Q  ,  and  period  of 
rise,  Tf .  The  resulting  equations  Ire  given  below: 

Qp  *  antilog  (1 .85*.6941ogA+,000099Pd)  (1) 

Tf  -  antilog  ( . 66* . 3591ogA- .000067Pd )  (2) 

where: 

Qp*peak  discharge  of  unit  hydrograph  in  CFS 

Tr*ti*e  of  rise  of  unit  hydrograph  in  hours 


A*drainage  area  in  square  mile?- 
Pd»populati°n  density  In  persons  per  sq.  mile 


The  error  in  using  equation  l  to  estimate  Q  lor 
the  13  basins  ranges  from  .  3X  to  33t,  with  an  average 
of  LOX.  The  correlation  coefficient  for  equation  1 
is  .991.  The  error  in  using  equation  2  to  estimate  T 
for  the  13  basins  ranges  from  42  to  952,  with  an 
average  of  242.  The  correlation  coefficient  for 
equation  2  is  .831. 

Equations  1  and  2  were  used  to  develop  unit 
hydrograph  parameters  for  the  ungaged  subha*.ins  in 
the  Red  Run  basin.  Required  population  densities  for 
the  years  1940,  1950,  I960,  and  1975  were  obtained 
from  census  data.  Population  projections  for  the 
years  2000  and  2025  were  obtained  fro*  the  South¬ 
eastern  Michigan  Council  of  Governments.  Tn  facili¬ 
tate  use  of  HEC-1 ,  values  for  Q  and  Tr  from  equations 
1  and  2  were  transformed  to  ^equivalent  values  of 
Snyder  unit  hydrograph  parameters,  t  and  C  . 
Adopted  values  for  these  parameters  art  shown  In 
Table  3. 
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Table  3.  Unit  Hydrograph  Parameters  and  Impervious  Areas  for  the  Red  Run  Drain 


Xsar 

1940 

1950 

I960 

1975 

2000 

r 

2025 

n 

m 

■ 

■ 

lap. 

■ 

B 

■ 

■ 

■ 

Basin 

Area 

mm 

■PH 

Ea 

C 

C 

lap. 

EH 

Hi 

Sq.  Hi. 

MM 

■flj 

Mi 

D 

B3 

P 

Hi 

P 

H 

13 

ID 

■3 

IWalva  Town 

37.13 

10.2 

.73 

.04 

7.3 

.83 

.07 

5.1 

.94 

.10 

4.6 

.96 

.10 

4.7 

.97 

.10 

4.6 

.97 

.11 

Henry  Crams 

21  .OS 

- 

- 

- 

- 

- 

- 

- 

- 

- 

8.2 

.70 

.05 

7.8 

.71 

.05 

7.3 

.71 

.05 

Local  Inflow 

at  Baar 

Crnak 

3.08 

7.7 

.43 

0 

7.5 

.52 

0 

5.0 

.61 

.04 

3.0 

.70 

.05 

2.3 

.75 

.09 

2.1 

.76 

.10 

Baar  Crack 

22.78 

13.7 

•58 

0 

13.5 

.58 

0 

9.1 

.66 

.04 

5.6 

.81 

.08 

5.1 

.84 

.08 

4.0 

.91 

.10 

Big  Baavar 

4.44 

14.3 

.38 

0 

14.1 

.59 

0 

13.6 

.59 

.01 

4.2 

.62 

.04 

2.5 

.71 

.08 

2.2 

.73 

.09 

Schoanharr 

12.83 

10.9 

.36 

0 

10.8 

.56 

0 

10.4 

.57 

.01 

4.7 

.76 

.09 

4.1 

.79 

.08 

3.5 

.82 

.09 

Local  Inflow 

at  Starling  - 
Balia  £ 

3.8 

.49 

0 

5.7 

.50 

0 

5.5 

.50 

.01 

3.6 

.60 

.05 

2.7 

.65 

.07 

2.2 

.68 

.08 

Starling  Dr. 

6.29 

- 

- 

- 

- 

- 

1 

- 

- 

- 

3.6 

.60 

.05 

3.9 

.68 

.06 

2.8 

.75 

.09 

Plua  Brook 

26.84 

16.0 

.60 

0 

15.7 

l _ 

.60 

B 

15.1 

.61 

.01 

12.7 

.60 

.03 

9.5 

.70 

.04 

9.1 

.71 

.04 

(1)  t  -  Snyder  "Standard"  Lag  in  Hour* 
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(2)  C  -  Snyder  Peak  Discharge  Coefficient 

P 

(3)  lap*  •  Hydro logically  significant  Impervious  area  expressed  ss 

a  proportion  of  subbasin  area 


Infiltration  Analysis  -  Hydrographs  recorded  for 
eight  subbasins  were  analyzed  to  determine  average 
soil  infiltration  capacity.  In  addition,  data  from 
other  subbasins  of  the  Clinton  River  and  the  adjacent 
River  Rouge  were  analyzed  and  published  by  Brater  and 
Sherrill  (2).  There  ia  considerable  scatter  around 
the  average  seasonal  values,  depending  primarily  on 
antecedent  rainfall  and  evapo-transpiration.  The 
adopted  average  values  for  infiltration  were  .40 
inch/hour  for  summer  conditions  and  .10  inch/hour  for 
winter  conditions. 

Impervious  Areas.  An  important  effect  of  urbaniza- 
tion  is  to  create  additional  impervious  area.  An 
estimate  for  the  area  of  impervious  surfaces  may  be 
made  from  aerial  photographs,  but  all  precipitation 
on  such  impervious  areas  does  not  contribute  to  sur¬ 
face  runoff.  For  this  reason,  Brater  (2)  employed 
the  concept  of  a  "Hydrological ly  Significant  Imper¬ 
vious  Area"  (HSIA),  which  represents  the  percent  of 
the  total  area  which  always  contributes  a  runoff 
equal  to  the  total  precipitation  minus  retention,  or, 
in  other  words,  almost  100Z  runoff.  Values  for  HSIA 
can  be  derived  from  the  analysis  of  hydrographs 
resulting  from  large  storms  occurring  after  a  long 
summer  dry  spell,  when  the  infiltration  capacity  of 
the  pervious  portion  is  very  high  and  therefore  the 
entire  surface  runoff  is  derived  from  the  Hydro- 
logically  Significant  Impervious  Area*  In  south¬ 
eastern  Michigan,  the  HSIA  is  approximately  twenty 
five  percent  of  the  gross  physical  impervious  area  as 
determined  from  aerial  photographs.  Brater  showed  a 
correlation  between  HSIA  and  population  density,  P. , 
for  twelve  watersheds  in  southeastern  Michigan  (2;. 
This  relationship  is  as  follows: 

HSIA  -  1.38  Pd  (3> 


where: 

=  the  population  density  in  thousands  of 

persons  per  square  mile 

HSIA  is  in  percent  of  total  area 

Prec ipi tat  ion  -  In  order  to  simulate  the  hydrologic 
response  of  the  Clinton  River  Basin  using  the  HEC-l 
program,  it  was  necessary  to  review  and  collect  the 
pertinent  precipitation  data  within  the  region  for  a 
number  of  historic  storms.  Precipitation  amounts 
were  input  into  the  HEC-l  program  on  a  subbasin  basis 
using  the  Ihiessen  polygon  or  isohyetal  method,  de¬ 
pending  on  the  extent  of  the  data  available. 

Stream  Flow  Routing  -  Stream  flow  routing  was  done 
using  the  Modified  Puls  methods.  With  the  Modified 
Puls  method,  discharge  is  considered  a  unique  func¬ 
tion  of  storage.  A  number  of  reaches  were  estab¬ 
lished  for  the  Main  Branch  of  the  Clinton  River  and 
the  Red  Run  Drain.  Storage  volumes  were  determined 
for  each  reach  using  tin?  HEC-2  Water  Surface  Profiles 
Program.  (It)  Discharge-storage  relationships  were 
then  developed  to  provide  the  required  input  for  the 
HEC-l  program. 

Calibration  and  Verification  -  Six  separate  IIKC-l 
data  models  were  developed  in  this  study,  r»*pr**- 
senting  different  stages  of  the  development  and  ur¬ 
banization  of  the  Clinton  Rivr*r  Watershed.  These 
node  Is  represent  conditions  in  the  1940's,  1950's, 

1940's,  1970's,  year  2000,  and  year  2025.  The  eight 
listoric  storm  events  listed  in  Tabic  4  were  simu¬ 
lated  to  verify  that  the  parameters  used  were 
reasonable.  From  gage  records,  flood  hydrographs 
were  determined  lor  thf*  Clinton  Rivr*r  at  Fraser, 
which  includes  runoff  from  both  the  Upper  Clinton 
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OISCHARGE  (CFS  *  102>  |  DISCHARGE  (CFS  x  «£) 


Figure  4.  Hydrogr;  Reproduction 


CLINTON  R.  AT  FRASER 
25  -  29  JUNE  1968 


Figure  5.  Hydrograph  Reproduction 


CLINTON  R.  AT  FRASER 
18  APRIL  -  22  APRIL  1975 


Figure  6.  Hydrograph  Reproduction 


River  and  the  Red  Run  Drain;  the  North  Branch  at  Mt. 
Clemens,  which  includes  runoff  from  the  North  Branch; 
and  the  Clinton  River  at  Mt.  Clemens,  which  includes 
runoff  from  almost  the  entire  Clinton  River  Basin. 
The  rainfall  associated  with  each  historic  event  was 
input  to  HEC-1  and  computed  flood  hydrographs  were 
compared  to  the  observed  hydrographs.  The  infiltra¬ 
tion  parameter  DLTKR,  which  is  an  index  of  antecedent 
moisture  conditions,  (8)  was  calibrated  within  rea¬ 
sonable  limits  to  provide  the  best  hydrograph  repro¬ 
duction.  Typical  hydrograph  reproductions  are  shown 
in  Figures  4,  5,  and  6. 


Frequency  Curves  At  Caged  Sites 

Two  key  stream  gage  sites  in  this  study  are  tb«> 
Clinton  River  near  Fraser  and  the  Clinton  River  at 
Mt.  Clemens.  The  Fraser  gage  has  been  in  operation 
since  1947  and  measures  runoff  from  a  444  square  mile 
area  that  includes  the  Upper  Clinton  River  and  the 
Red  Run  Drain.  The  Mt .  Clemens  gage  has  been  in 
operation  since  1934  and  measures  runoff  from  714 
square  miles,  which  is  nearly  the  entire  Clinton 
River  basin.  Records  from  these  two  gages  were  of 
substantial  importance  in  this  study  for  calibrating 
the  rainfall-runoff  simulation  model ,  and  for  assign¬ 
ing  the  exceedance  frequencies  to  synthetic  storm- 
loss  rate  combinations. 
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Table  4. 

Historic  Storm 
Verification  an 

Events  f) soil  lor 
if  Calibration 

Observed 

Obso  rved 

Peak  Discharge 

Peak  Discharge 

Date  of 

at  Fraser 

at  Mt.  Clemens 

Storm 

(cfs) 

(cfs) 

29  August  - 
2  September 
1975 

4,400 

5,200 

19-22  April 
1975 

5,500 

12,100 

4-9  March 

1974 

3,800 

10,200 

26-30  June 

1973 

4,400 

3,900 

25-29  June 

1968 

7,000 

1  1,000 

27  April  - 
2  May  1956 

5,200 

13,500 

10-13  May  1948 

8,000 

14,500 

3-8  April  1947 

- 

21,200 

IIm*  ^nes  <*1  annual  peak  d  i  scti/jrg<*s  ai 

pacli  site  is  nuns  tat  i  unary  and  nonhomugeneous  -  noil- 
stationary  because  of  the  time-dependent  effects  of 
urban  development;  nonhomogeneous  because  in  some 
years  a  winter  event  produces  the  annual  maximum 
discharge,  whereas  in  other  years  a  summer  event 
produces  the  annual  maximum.  Winter  runoff  events 
are  characterized  by  cyclonic  storms  and  very  low 
loss  r.i»<*s;  siimmir  events  by  t  bunder  si  u  rms  and 
relatively  large  loss  rates.  The  following  pro¬ 
cedures  were  used  for  I  runs  /  ormi  ng  tin  recorded 
annual  peak  discharges  at  Kras*  r  and  Mi  .  Clemens 
into  i  tins  i  stent  series  that  ref  let  ts  existing  land 
use  condition^  and  channel  cha  r.ii  ter  i  t  i»s  (  i  2  )  . 

As  indicated  previously,  input  data  sets  for 
HKC-1  were  developed  and  calibrated  to  reflect  exist¬ 
ing  watershed  conditions  and  conditions  in  1960, 
1950,  and  1940.  Synthetic  winter  and  summer  storm 
events  were  prepared  using  storm  patterns  (i.e.,  time 
distributions)  developed  by  Brater  and  Sherrill  (2). 
They  found  that  *. 

"Typical  large  rainstorms  in  Southeastern 
Michigan  have  a  nearly  symmetrical  time 
distribution  with  the  maximum  hour  near  the 
center  of  the  rain  storm.  In  winter  the 
maximum  hour  contains  about  24  percent  of 
the  rain  and  in  summer  the  maximum  hour 
usually  provides  about  55  percent  of  the 
total  storm  rainfall  .  . 

The  magnitude  of  a  "base"  storm  was  fixed  by 
adopting  a  25-year  recurrence  interval,  24-hour  dur- 


Figure  7.  Teak  Discharges  v»  Rat  ion  of  Kahtl.dl  at  the  i*  i  »a  i.ng*' 
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at  ion  piecipitat  ion  depth  from  NOAA's  Technical  Paper 
Ho.  40  (9).  The  recurrence  interval  of  25  years  was 
arbitrary.  It  was  not  assumed  in  the  approach  des¬ 
cribed  herein  that  the  frequency  of  runoff  is  equal 
to  the  frequency  associated  with  the  synthetic  storm. 
Rather,  the  purpose  was  to  establish  a  base  storm  for 
use  in  subsequent  steps  in  the  analysis. 

Six  synthetic  storms  for  each  seaaon,  summer  and 
winter,  were  developed  by  applying  scaling  ratios  to 
the  base  storms.  Ratios  were  selected  so  that  peak 
discharges  resulting  from  the  storms  would  cover  the 
range  of  discharges  required  to  define  discharge- 
frequency  curves  at  the  Fraser  and  Mt .  Clemens  gage 
locations.  Winter  synthetic  storms  were  developed 
using  scaling  ratios  of  .2,  .6,  1.0,  1.4,  1.8  and  2.2. 
Suaer  synthetic  storms  were  developed  with  ratios  of 
.25,  .5,  1.0,  1.5,  2.0,  and  2.5. 

The  mul ti plan-mu  1 1 i rat io  capability  of  HEC-1 
permitted  calculation  of  discharge  hydrographs  for 
six  storms  and  four  sets  of  runoff  parameters  (for 
watershed  conditions  in  1975,  I960,  1950,  and  1940) 
in  one  run  of  the  program.  Hence  two  runs  were 
required  -  one  for  summer  conditions  and  one  for 
winter  conditions.  Resulting  relationships  for  cal¬ 
culated  peak  discharge  versus  scaling  ratio  for  the 
Fraser  and  Mt.  Clemens  gages  are  shown  in  Figures  7 
and  8.  The  relationships  were  used  to  transform  each 
recorded  annual  peak  discharge  (for  each  season)  to 
an  exist ing-cond it  ion  annual  peak.  This  was  achieved 
by  entering  the  appropriate  set  of  curves  horizon¬ 
tally  with  the  recorded  peak  discharge  to  the  inter¬ 
section  with  a  curve  (or  interpolated  curve)  for  the 


year  in  which  the  discharge  was  recorded,  and  then 
moving  vertically  to  an  intersection  with  the  1975 
curve,  from  which  the  transformed  peak  discharge 
could  be  determined.  The  largpBt  of  the  two  trans¬ 
formed  peaks  for  each  year  was  incorporated  in  a  tint- 
series  of  annual  peaks.  Statistical  analyses  of  the 
resulting  series  was  performed  using  techniques  in 
Bulletin  17  of  the  Water  Resources  Council  (10).  The 
resulting  frequency  curves  are  shown  in  Figures  9  and 
10.  Also  shown  are  frequency  curves  based  on  ob¬ 
served  (untrans formed)  peak  discharges.  Urbanization 
has  had  a  substantial  influence  on  the  frequency 
curve  for  Fraser,  whereas  the  influence  at  Mt . 
Clemens  is  minor. 


Seasonal  Frequency  Curves  At  Fraser 

Subsequent  development  of  discharge  frequency 
curves  for  ungaged  sites  required  separate  frequency 
curves  for  winter  and  summer  conditions  for  the 
Clinton  River  at  Fraser.  The  curves  are  shown  in 
Figure  11.  The  annual  curve  is  the  same  as  shown  in 
Figure  9.  The  winter  curve  was  obtained  graphically 
using  Weibutl  plotting  positions  for  the  transformed 
winter  peaks.  The  series  of  summer  peaks  was  trunca¬ 
ted  because  for  a  number  of  years  the  observed  sunnier 
peaks  were  smaller  than  the  base  value  used  by  the 
U.S.  Geological  Survey  aa  a  basis  for  publishing 
instantaneous  peaks.  The  summer  frequency  curve  was 
therefore  obtained  by  subtraction  using  the  following 
equation: 


Figure  8.  Peak  Discharge;*  vs  Ratios  of  Tula  I  Rainfall  at  the  Mi.  Clemons  Ongv 
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Figure  9.  Frequency  Curve  at  Fraser 


PANNUAL  -  PWINIER  (4) 

PSUMMER  -  - 

1  “  PWINTBR 

where: 

P_lfl___0  *  Exceedance  probability  for  a  given 

D  Urlrul  K  •  <  *  .  -• 

discharge  on  the  summer-season  fre¬ 
quency  curve. 

PWINTER  *  Exceedance  probability  for  a  given 

discharge  on  the  winter-season  fre¬ 
quency  curve. 

PANNUAL  *  Exceedance  probability  for  a  given 

discharge  on  the  annual  (all-season) 
frequency  curve. 


The  available  transformed  summer  peaks  are  plotted  in 
Figure  11  for  comparison  using  Weibull  plotting 
posit  ions. 

Frequency  Curves  At  Ungiged  Sites 

The  basis  for  developing  peak-d ischarge  fre¬ 
quency  curves  at  ungaged  sites  throughout  the  140 
square  wile  Red  Run  Basin  was  to  use  balanced 
synthetic  storms  as  inputs  to  calibrated  HEC-l  models 
for  the  basin*  Synthetic  storm-loss  rate  combina¬ 
tions  were  assigned  exceedance  frequencies  in  ac¬ 
cordance  with  exceedance  frequencies  of  calculated 
peak  discharges  obtained  from  the  previously  develop¬ 
ed  frequency  curves  for  Fraser. 


Synthetic  winter  and  summer  storms  were  based  on 
criteria  from  NOAA's  Technical  Paper  No.  40  (9)  and 
the  time  distributions  developed  by  Brater  and 
Sherrill  (2).  Application  of  areal  adjustments  from 
Technical  Paper  No.  40  to  point  rainfalls  of  various 
durations  results  in  storm  hyetographs  that  are  a 
function  of  drainage  area  size.  Therefore,  separate 
HEC-l  runs  were  made  using  storm  hyetographs  as¬ 
sociated  with  drainage  areas  of  10  square  miles  and 
400  square  miles.  Peak  discharge  for  any  particular 
(ungaged)  drainage  area  size  was  obtained  by  log¬ 
arithmic  interpolation  between  the  peak  discharge 
associated  with  the  storm  hyetographs  for  10  square 
miles  and  400  square  miles. 

Winter  and  summer  frequency  curves  were  develop¬ 
ed  at  all  ungaged  locations  of  interest,  for  existing 
land  use  conditions  and  for  land  use  conditions  in 
the  years  2000  and  2025,  using  the  techniques  des¬ 
cribed  above.  Winter  and  summer  frequency  curves  for 
a  given  land  use  condition  were  then  combined  to 
obtain  .in  annual  curve  based  on  the  probability  of 
union  of  independent  events.  Figure  12  shows  a 
typical  set  of  Frequency  curves  foe  an  engaged  loca* 
t  ion. 


Kindi ngs 

The  Rod  Run  Drain  and  bower  Clinton  River  water¬ 
sheds  have  undergone  urbanization  since  th *  1940's. 
The  effect  of  thi  s  urbanization  on  the  peak  dis¬ 
charges  recorded  at  the  1  user  and  Mt .  Clemens  gages 
can  be  seen  in  Figures  7  and  ft.  Urbanization  has  had 
a  substantial  influence  on  the  frequency  »»{  d isobars*- 
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for  Fraser,  as  shown  in  Figure  9.  The  peak  discharge 
at  Fraser  is  comprised  almost  entirely  of  runoff  from 
the  Red  Run  Drain,  which  peaks  much  sooner  than 
runoff  from  the  main  branch  of  the  Clinton  River.  The 
Red  Run  has  been  substantially  urbanized,  and  so  this 
change  in  frequency  ia  not  unexpected.  The  influence 
of  urbanisation  at  Ht.  Clement  is  minor  as  shown  in 
Figures  8  and  10.  Peak  discharges  at  Mt.  Clemens  are 
a  combination  of  runoff  from  the  North  Branch  Water¬ 
shed,  largely  rural,  and  the  Red  Run  Drain.  Urbani¬ 
sation  has  decreased  the  peaking  time  of  the  Red  Run 
Drain,  which  tends  to  move  apart  the  time  of  occur¬ 
rence  of  the  peek  runoff  rates  of  the  North  Branch  and 
the  Red  Run  Drain.  This  phenomenon  has  tended  to 
counteract  the  increase  in  peak  discharge  rates  due 
to  urbanisation. 

Summary 

The  development  of  peak-discharge  frequency 
estimates  at  ungaged  locations  in  an  urbanising 
watershed  for  existing  and  future  conditions  is,  in 
our  opinion,  one  of  Che  most  challenging  tasks  a 
hydrologist  can  face.  Three  factors  aided  substan¬ 
tially  in  the  conduct  of  the  study;  (1)  availability 
of  results  of  previous  extensive  studies  pertaining 
to  the  hydrology  of  Southeastern  Michigan,  (2)  a 
relative  abundance  of  precipitation  data  and  O)  the 
finding  of  a  statistically  significant  index  of  ur¬ 
banization,  population  density,  to  which  runoff  para¬ 
meters  could  be  related  in  a  convenient  fashion. 


The  basic  approach  was  to  use  a  single  event 
rainfali-runof f  model  to  transform  nonstationary, 
nonhoisogencous  series  of  recorded  annual  peak  dis¬ 
charges  to  consistent  series  representing  existing 
and  future  conditions.  The  problem  of  nonhotaogcneity 
was  handled  by  using  separate  analyses  for  the  summer 
and  winter  seasons.  Norist  at  ionari  ty  was  treated  by 
developing  a  set  of  storm  magnitude  vs.  peak  dis¬ 
charge  relationships  for  vsrious  points  in  time  and 
using  these  relationships  to  transform  nous tat  ionary 
series  to  stationary  ones.  In  utilising  synthetic 
storms  to  develop  peak  discharges  for  ungaged  loca¬ 
tions,  storm-loss  rate  combinations  were  assigned 
exceedance  frequencies  equal  to  the  exceedance  fre¬ 
quencies  of  calculated  peak  discharges  at  the  gage  at 
Fraser. 
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Figure  10.  Frequency  Curve  at  Mt.  Clemens 
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Figure  11.  Seasonal  Frequency  Curves  at  Fraser 


Figure  12.  Typical  Kronuetii  y  Curves  .11  .m  l«*» 
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